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We investigate the behavior of the critical charge for spontaneous pair production, Zc, defined 
as the charge at which the total energy of a A'-shell electron is i? = —rUe, as a function of the 
radius R of the charge distribution. Our approach is to solve the Dirac equation for a potential 
V{r) consisting of a spherically symmetrical charge distribution of radius R and a Coulomb tail. For 
a spherical shell distribution of the type usually associated with color-flavor locked strange quark 
nuggets, we conflrm the relation Zc = 0.71 i?(fm) for sufficiently large R obtained by Madsen, 
who used an approach based on the Thomas-Fermi model. We also present results for a uniformly 
charged sphere and again find that Zc ^ R for large enough R. Also discussed is the behavior of 
Zc when simple ad hoc modifications are made to the potential for < r < _R. 

PACS numbers: 21.65. Qr,24.85.-|-p 



1. INTRODUCTION 

In a recent paper we examined the changes that develop on strange quark nuggets (SQN's) in space 
environments and discussed various signatures based on the charges they would acquire in these circum- 
stances. Among these was the possibility that the surface charge on the nugget might exceed the critical 
charge Zc associated with the onset of spontaneous pair production. The condition for this to occur, that 
the energy of a K-sheW electron satisfies E = — TOc, depends on the radius of the charge distribution 0, S3- 
For color-flavor locked SQN's it has been shown that the charge and the baryon number are related as 
Zq = 0.3 A^/^ ~ i?^ and the question arises as to whether the nugget charge Z^ ever exceeds the critical 
charge for any value of R. This issue was examined in a paper by Madsen 0], who solved the Thomas- 
Fermi model for the case of the critical energy E = —me. He showed that the color-flavor locked relation 
Zc ~ Zq = 0.3 A^/^ is valid for nugget radii up to about 40 fni at which point there is a transition to a 
relation between the critical charge Zc and the radius R given by Zca = 2meR or Zc = 0.71 i?(fm) for 
R > 400 fm. 

Our approach to the determination of the critical charge for a spherically symmetrical distribution is one 
of solving the Dirac equation with the energy fixed at the critical value E = —m^. By matching the ratios 
of the radial solutions for r < R and R > r at the boundary, it is then possible to determine Zc for a given 
R to any desired accuracy. In this way, the variation of Zc as a function of R can be obtained. For the 
spherical shell, we find that the relation Zcct = 2meR is satisfied for R > 800 fm. For R < 200 fm, Zc is 
slowly varying with respect to R. In Section [51 we derive the dependence of Zc on R for two spherically 
symmetrical distributions, a charged shell and a uniformly charged sphere. This is followed, in Section [31 by 
a discussion and some conclusions. 
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2. EVALUATING THE CRITICAL CHARGE FOR SPHERICAL DISTRIBUTIONS 

In the following cases, we treat potentials of the form 

r V(r) for 0<r <R 
[ for r > R 

where V(r) takes on different forms depending on the particular charge distribution. To determine the critical 
charge, we solve the Dirac equation, whose upper and lower components G(r) and F{r) satisfy Q 

= --G{r) + {m + E-V{r))F{r) (2a) 
dr r 

^ -F{r) + {m-E + V{r))G{r), (2b) 
ar r 

in the special case E = — m. Here k = £ for j = ^ — i and k = —{£ + 1) for j ^ {i + |). 

Depending on V(r), the form of the interior solution {0 < r < R) varies and, in the case of the uniformly 
charged sphere, must be obtained numerically. The form of the exterior (r > R) solution is known [ii]. By 
setting E = —m in Eqs. ([2]) and eliminating F{r) in favor of G(r), one obtains the differential equation 

x^^G{x) + x-^G{x) - x''G{x) + 4(^2^2 _ k')G{x) = , (3) 
dx'^ dx 



with x — \/8Zarm. The solution to Eq. ([3]) is a modified Besscl fmiction of imaginary order 



G{x)^Ki^{x) with u = 2y/Z^a^ - . (4) 

To determine the critical charge, we use the condition that the ratio F{r)/G{r) for the exterior and interior 
solutions must be equal at the boundary r = i? For the exterior solution, this ratio can be obtained 
using Eq. (|4]) and Eq. ((2a|) with E = — m. The result is 



(5) 



F(r) _ K r dG{r) _ k VSZamr KUVsZ amr) 

GM~^ ZaG{r) dr " 2Za K,^{y/SZamr) ' 

2.1. Spherical Shell 

For a constant interior potential, V(r) = Vq, G{r) satisfies 

G"{r) + (Vo{2m + K) - G(r) = . (6) 

In particular, the interior potential for a uniformly charged spherical shell of radius R is 

V(.) = -f, (7) 

and, taking x — 2mr, Eq. ^ becomes 

G"{x)+(k'-^-^^)G{x). (8) 
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Here, we have used the fact that, for interior solutions which remain finite at r — 0, k{k+1) — i{£+l). The 
constant is 



i?2 Vi?2 



(9) 



and R = y/2mR. The solution to Eq. ([8]) is G{x) = xji{kx). Again, using Eq. ((2a)l . the ratio of the interior 
solutions is 



G{r) Zar Zar 



je{2mkr) 



(10) 



The critical value of Za can now be determined by requiring that Eqs. ([5]) and (jlOp are equal at r = R. The 
resulting condition is 



a;(x- i?2) 
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itV2; — ^ 



(11) 



where x = Za. Accurate numerical evaluations of the Bessel functions in Eq. (jlip are available in both the 
Mathematica and Maple packages. The location of the critical value of Za can be estimated by plotting 
both sides of Eq. (fTTj) and then using a search routine to obtain the precise value. This is illustrated in FiglT] 
for i? = 12fm (i? = 0.25) and £ = (k = —1) which corresponds to the s^j^ states. The lowest crossing, 
Za = 1.296 in this case determines the critical charge. As R increases, the first solution continues to 





FIG. 1: The solutions for the critical value of Za corre- 
sponding to the X-shoU is shown. The red curves are the 
exterior solution for R = 0.25 and the green curves are 
the corresponding interior solution. 



FIG. 2: The same as Fig[T]for R = 6.0. 



increase, approaching the value Za = i?^ — 2mR. This is illustrated in Fig. [2] for R ~ 6 and shown in detail 
in Table H 

If R is expressed in femtometers, the relation between the critical charge Zq and R can be written as 



Zc = 0.71 R (fm) , 



(12) 



which is the result found by Madsen 0], who derived it in a totally different way (Thomas-Fermi screening). 
He argues that Eq. (11211 holds for R > 400 fm whereas we find that onset of validity is more like R > 
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£ = 


£ = 1 


£=1 




K = -1 


K = 1 


K = —2 


R 








0.10 


1.154 


1.207 


2.120 


0.25 


1.296 


1.414 


2.251 


0.50 


1.557 


1.792 


2.505 


1.00 


2.275 


2.730 


3.206 


2.00 


4.925 


5.560 


5.717 


3.00 


9.608 


10.147 


10.191 


4.00 


16.411 


16.811 


16.824 


5.00 


25.291 


25.586 


25.590 


6.00 


36.216 


36.437 


36.439 


7.00 


49.166 


49.337 


49.336 


8.00 


64.131 


64.267 


64.267 


9.00 


81.106 


81.216 


81.216 


10.00 


100.087 


100.178 


100.178 



TABLE I: The dependence of Zcoi on R = V2mR 
is shown for the spherical shell distribution. The 
columns correspond to the 5^/2 {^^ = Pl/2 
{k = 1) and P3/2 = —2) states. 
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FIG. 3: The critical value of Z for the spherical shell obtained 
by solving the Dirac equation (solid green line), the approach 
to Eq. I II2I 1 (red dashed line) and the expression for charge of 
the color-flavor locked SQN, Zq = O.SA^/S = 0.726 (blue 
dot-dashed line) are shown. 



800 — 1000 fm, as can be seen in Fig.[3l In this figure, we show tlic critical value of Z obtained by solving the 
Dirac equation for k = — 1 (solid green line) and the approach to Eq. (red dashed line) together with 
the expression for charge of the color-flavor locked SQN, Zq = 0.3 A^/^ = 0.726 i?^ (blue dashed hue). 

For SNQ's with R > 20 fm, the value of Z obtained from the relation Zq = 0.3^^/"^ exceeds the critical 
value and the nugget will screen its charge by producing pairs until its charge falls to the green curve in 
Fig.m As the nugget radius increases, the screening process continues to limit the nugget's charge to the 
values on the critical line. Ultimately, the relation between the critical charge and the nugget radius is given 
by Eq. p2)) . Beyond some radius, Madsen argues that the time required to screen the charge will be such 
that the relation between the charge and the radius will again be of the form Z ^ . The analysis leading 
to this conclusion is quite unrelated to the task of determining the critical charge using the Dirac equation 
or the Thomas- Fermi model, and we have nothing to add in this regard beyond what is discussed in Ref. [l| . 



2.2. Uniformly Charged Sphere 



For a uniformly charged sphere of radius i?, the interior potential is 

Za f r"^ 



(13) 



Inserting this potential into Eqs. ([2|) with E = — m and letting y = r/ R, the functions G{y) and F{y) satisfy 



dG{y) 

dy 
dFjy) 

dy 



y ^ 

'-F{y) + (R' + l{y^-Z))G{y)^ 



(14a) 
(14b) 
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where, as before, R ~ \/2inR and x ~ Za. In this case, the ratio of F{r)/G{r) at the boundary r ^ R 
{y = 1) must be obtained by numericaUy integrating Eqs. (jl4[) for a given R and a range of Za. This can 
be done using routines available in Mathematica or Maple and the result compared with the exterior ratio 
Eq. ^ to obtain an estimate of the lowest value of Za. The value can be refined to any desired accuracy 
by searching in the vicinity of the estimated solution. 
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K = -2 


R 


Zcce 
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Zcce 


0.10 
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3.019 


2.00 
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3.00 
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54.00 


55.95 


10.00 


67.83 
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12.00 
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96.00 
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14.00 


131.83 


130.67 


132.60 


16.00 


171.83 


170.67 


172.60 


18.00 


217.16 


216.00 


221.02 


20.00 


267.82 


266.67 


268.60 



T ABL E II: The dependence oi Zcce on R = 
\/2mR is shown for the uniform charge distribu- 
tion. The columns correspond to the Sxji {k = 
— 1), pi/2 (t = 1) and P3/2 (k = —2) states. 
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FIG. 4: The critical Z for the uniformly charged sphere of 
radius K is shown in the solid green line. The dashed red line is 
Z = 0.475 R. 



Values of Zca for a range of R are given in Table HIl and the critical charge for k = — 1 is shown as the solid 
green line in Fig. 3) The dependence of the critical charge on the radius is quite similar to the spherical shell 
case shown in Fig. [3] in that there is a region of slow variation followed by a transition to a linear dependence 
on R. Overall, the values of Zq for the uniformly charge sphere are lower at a given R than in the spherical 
shell case, and the asymptotic i?-dependence is Zca ~ ^m.f.R. 



3. DISCUSSION AND CONCLUSION 



We have obtained the critical charge as a function of the radius for two spherically symmetrical charge 
distributions by solving the Dirac equation to determine the values of Za that produce /^-shell and L-shell 
electron energies of i? = — m. In both cases, the spherical shell and the uniformly charged sphere, the 
dependence of the critical charge on distribution's radius is similar. 

For the shell distribution, we exactly reproduce the large R relation Zca = 2mR found by Madsen Q, 
who used an approach based on the Thomas-Fermi atomic model. We see this behavior for s-states as well 
as for p-states with j ~ 1/2 and 3/2. For radii of order a few fcrmis. the behavior of the critical value of Za 
is controlled by the order of the modified Bessel function in Eq. v = Z'^cP- — k?- . Unless Za > \k\, 
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the Bcsscl function is not oscillatory and there is no solution to the boundary condition. Hence, the critical 
Zo. begins at a value slightly larger than |k| for small i?, and evolves to the line Za = 2mR as R increases. 
This behavior can be seen in TableU and is illustrated in Fig. ([5]) . 
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FIG. 5: The behavior of the critical values of Z are shown as a function of the shell distribution radius R for levels with k = — 1 
and K = —2 (solid green lines) as they approach the line Za = 2mR (red dashed line). The SQN charge, Zq = 0.3 A^/^, is also 
shown (blue dot-dashed line). 

The situation for the uniformly charged sphere is quite similar. Generally speaking, the critical values of 
Za for a given radius are slightly smaller than the corresponding values for the charged shell. For large R, 
the critical value of Za is again proportional to R with a smaller slope, |to. The details are given in Table 

IE 

In the discussion thus far. we examined the i?-dependence of the lowest critical value of Za for a given 
value of K and found a universal linear dependence for sufficiently large values of R. This behavior persists 
when higher critical values of Za for a particular k are examined, as illustrated in FiglHl where we plot the 
evolution of the first and second critical values of Z for k = — 1 as a function of R in the spherical shell case. 
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FIG. 6: The behavior of the first and second critical values of Z for k = — 1 are shown as a function of the spherical shell 
radius R (solid green lines) as they approach the line Za = 2mR (red dashed line). The SQN charge, Zq = 0.3 A^/^, is also 
shown (blue dot-dashed line). 

In fact, all states will have a Zca dependence on R with the same general shape as the curves in FiglBJ they 
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will simply start at larger values on the vertical (Z) axis. We used Eq. at i? = 2fm with k from —1 to 
3 {£ from to 3) to do a count and found 64 states with Zca < 4. These range from the lsi/2 state with 
Zc = 1-15 to the 11/5/2 state with Zca = 3.97 and can accommodate 472 electrons. 

We also looked at an ad hoc modification of the interior potential , 

that is less negative at r = by a factor of two than the shell potential of Eq. ([7]), rather than being more 
negative as is the case for the uniform sphere potential of Eq. The resulting values of Zccx are slightly 
larger than those in Table H] and have a very similar linear behavior for large enough R. 

Zc is the value of Zpf such that the least tightly bound electron in the ground state of the system has a 
binding energy of —2me- If that electron is removed, it becomes energetically favorable for the vacuum to 
produce an e^e~ pair. The interaction of the least bound electron with the other electrons is not taken into 
account in our calculations, but this should not change things by a great deal. 

To summarize, solving the Dirac equation with a spherically symmetrical charge distribution of total charge 
Ze and radius R for the energy E = —m predicts that the critical charge obeys Zc ^ R for sufficiently 
large R. This result holds for any value of the parameter k as well as for any of the critical values of Za 
associated with a given k. However, the approach to the linear dependence on R and its slope both depend 
on the particular charge distribution. 
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